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ABSTRACT 



The Brans-Dicke model with a variable cosmological term (BDA) has been investigated with use of the coupling constant of toi = lO'*. 
Parameters inherent in this model are constrained from comparison between Big Bang nucleosynthesis and the observed abundances. 
Furthermore, the magnitude redshift (m - z) relations are studied for BDA with and without another constant cosmological term in 
a flat universe. Observational data of Type la Supernovae are used in the redshift range of 0.01 < z < 2. It is found that our model 
with energy density of the constant cosmological term with the value of 0.7 can explain the SNIa observations, though the model 
parameters are insensitive to the m- z relation. 
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1. Introduction 

Astronomical observations indicate that the cosmological con- 
stant in the very early universe exceeds the present value by 
some 120 orders of magnitude, which is estimated in modern 
theories of elementary particles (Weinberg 1989). This is one of 
the fine tuning problems in cosmology called by the "cosmolog- 
, ical constant problem". To explain the puzzle in cosmology, new 
■ modified theories are needed beyond the standard model. That 
' behavior of the cosmological term has motivated various func- 
tional forms to the cosmological term. The mechanism to the 
dynamical reduction of the cosmological term is formulated as a 
time dependent function (Silviera & Waga ,1997j and in terms of 
a scalar field (Weinberg [19891 Huterer & Turner [T999] l. On the 
other hand, generalized scalar tensor theories have been investi- 
gated (Wagoner 1970, Endo & Fukui 1977 , Fukui et al. 2001 ). 
Among them a Brans-Dicke (BD) theory with a variable cos- 
' mological term (A) as a function of scalar field ((p) (Endo & 
Fukui 1 19771 1 has been proposed. This model has been inves- 
tigated for the early universe and constrained from Big Bang 
nucleosynthesis (BBN) (Arai et al. .1987. Etoh et al. [T9971 
Nakamura et al. I2006I I for the coupUng constant ai < 500. 
Present observations suggest that the value of w exceeds 40,000 
(Berti et al. 2003 Bertotti et al. 2005). Therefore, it is worth- 
while to reconstrain the parameters in the Brans-Dicke model 
with a variable cosmological term (BDA) for a new value of 
CO. BDA has played a very important role to explain the char- 
acteristics of the early universe (Arai et al. 119871 Etoh et al. 
[1997 Nakamura et al. I2006I I. However, an answer is needed 
to the question "How this model works at the present epoch?". 
Therefore we adopt the magnitude redshift (m - z) relations of 



Type la Supernova (SNIa) observations. This is because, the cos- 
mological term affects the cosmic expansion rate of the universe 
significantly even at the low redshifts. SNIa observations im- 
ply that the universe is accelerating around the present epoch 
(Perlmutter et al. 1999 Riess et al. |1998l |1999). 

In Sect. 2 formulation of BDA is reviewed. Parameters in- 
herent in this model are constrained in Sect. 3 from Big Bang 
nucleosynthesis for w - 10"^. In Sect. 4 the m - z relation is 
investigated for BDA with including another constant cosmo- 
logical term in a flat universe. Recent SNIa observational data 
(Astier et al. .2006. Riess et al. [20071 Kessler et al. [2009] l are 
adopted to constrain the models. Concluding remarks are given 
in Sect. 5. 



2. Brans-Dicke model with a variable cosmological 
term 

The field equations for BDA are written as follows (Arai et al. 
1 87r w / 1 „' 
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where is the scalar field and T^y is the energy-momentum ten- 
sor of the matter field. The Robertson- Walker metric for homo- 
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geneous and isotropic universe is written as (Weinberg II 9721) : 



ds^ 



-df + a (ty 



l-kr^ 



+ r^d9^ + r^sir?ed^'^ 



(3) 



where a{t) is the scale factor and k is the curvature constant. Here 
we adopt c - I. The expansion is described by the following 
equation derived from the (0, 0) component of Eq. ([T]i: 



laV onp k A cl) (p\ a (/> 
where p is the energy density. 
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We adopt the simplest case of the coupUng between the 
scalar and matter field is 
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where // is a constant. Assuming a perfect fluid for T^v, Eq. (|5]l 
reduces to the following: 
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where p is the pressure. 
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Fig. 1. Time evolutions of the scale factor in BDA with fi - 0.6 



A particular solution of Eq. ^ is obtained from Eqs. 0]) and which are compared to the Friedmann model. 
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(7) 



where p„,„ is the matter density at the present epoch. 

The gravitational "constant" G is expressed as follows 
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The density p and the pressure p are replaced by 

P - Pm+ Py, 

P ^ Py^ Pyl'i 

where the energy density of matter varies as p^ = Pmoa 



(9) 
(10) 

I The 
energy density of radiation is written as py — Py„a^'^ except e=^ 
epoch: py - prad + Pv + Pf± at f < Is, where subscripts rad, 
V and e=^ are for photons, neutrinos and electron-positrons, 
respectively (Nakamura et al. l2006l l. Subscript "0" indicates the 
values at the present epoch. 



Then, Eq. ^ is integrated to give 
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where B is an integral constant and here we use the normalized 
value of B: B* = B/(10"24g ^ ^,^-T>y 

Original Brans-Dicke theory is deduced for p - I and it is 
reduced to the Friedmann model when (p - constant and w » 1 . 
Physical parameters have been used to solve Eqs. (JUi, (|7]i and 
dnii: Go = 6. 6726 X IQ-^^cm^g-'s-^ //,) = 71 km s' ' Mpc^i 
(Spe rgel et al. l2003l l. and u = 10'' (Berti et al. l2003l Bertotti 
et al. 12005 1). Figure [T| shows the evolution of the scale factor in 
BDA for the several values of B*. We identify considerable devi- 
ations in BDA from the Friedmann model at f < 100 s, which de- 
pends on the specific parameters. Therefore BDA should be con- 
strai ned fro m BBN (Aral et al. [T987l Etoh et al. [T997l Nakamura 
etal.|2006l). 



3. Parameters constrained from Big Bang 
nucleosynthesis 

Big Bang nucleosyntheis provides powerful constraints on 
possible deviation from the standard cosmology (Malaney & 
Mathews. [1993]). As shown in Fig.lT] the expansion rates of BDA 
differs significantly from that of the standard Friedmann model. 

The abundance of light elements in BDA has akeady been 
investigated (Aral et al. 1987, Etoh et al. II 9971 Nakamura et al. 
2006). In the previous studies, the parameters inherent in BDA 
have been constrained for oj = 500. But we consider the case 
w = 10"* for convenience, because the Cassini measurements of 
the Shapiro time delay indicate w > 4 x 10"* (Berti et al. 2003) 
Bertotti et al. 2005 ). The detailed method of nucleosynthesis is 
described in Nakamura et al. 120061 

Figure |2] shows the calculated abundances of '*He, D, 
and ^Li for B* - 2 and p - 0.6. The +2cr uncertainties 
in nuclear reaction rates are indicated by the dashed lines. 
The horizontal dotted lines indicate the observational values 
of '♦He, D/H, and ^Li/Has follows: F^ = 0.2516 ± 0.0080 
(Fukugita & Kawasaki |2006|) . Y,, = 0.326 + 0.075 (Komatsu 
et al. l20T0b . D/H = (2.82 ± 0.21) x 10"^ (Pettini et al. 12008] ), 
^Li/H = (2.34 ± 0.32) X 10-1° (Melendez & Ramirez. l2004l ). 
Here two observational values of ''He are used. The solid verti- 
cal lines indicates the WMAP constraint of the baryon-to-photon 
ratio, // = (6.19 + 0.15) x 10""' (Komatsu et al. l2UT0l ). 

The intersection range of the two observational values of ''He 
is used to constrain the parameters. It is found that the values 
of 77 derived from ''He and D/H are tightly consistent with the 
value by WMAP, though the lower limit of ^Li/H is barely con- 
sistent. These agreements lead us to obtain the parameter ranges 
of 0.0 < yu < 0.6 and -2 < B* < 2. 
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Fig. 2. Light element abundances of ^He, D, and 'Li vs. 77 for 
B* - 2,fi - 0.6, and co - lO'*. Dashed lines indicate the +2cr un- 
certainties in nuclear reaction rates in each abundance. The hor- 
izontal dotted lines indicate the regions of observational abun- 
dances. The solid vertical lines indicate the baryon-to photon 
ratio rj. 

4. m- z relation in BDA with and without a constant 
cosmological term 

The distance modulus ju,/, of the source at the redshift z is 

//„, ^m-M^ 51ogio [(1 + z) n] + 25, (12) 

where m and M, are the apparent and absolute magnitudes, re- 
spectively and ri stands for the radial distance in units of Mpc. 

We adopt the SNIa (Astier et al. l2006l Riess et al. I^OTl 
Kessler et al. 120091 1 for which ;^f^ is defined by 



2 Y^ (MthJ - f^obsj) 
' =>' ^^ 






(13) 



where ju,/, , is given by Eq. ( fT2] ). ^ohs,i and o", are the observed 
values of distance modulus and their uncertinities. 

For the homogeneous and isotropic universe, the relation be- 
tween the radial distance and the redshift is derived from the 
Robertson-Walker metric as (Weinberg 120081) 



^- , f r'^'^sin^M V^r,) k^+\, 

r * ; n 
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where H - aja is the expansion rate written from Eq. Q as 



// = 






-2I '''' 



We conclude from the WMAP results that we live in a closely 
geometrically flat universe (Dunckley et al. 120091) . The present 
matter density p,„^ is obtained from Eq. Q as 



«-5(^-» 



-■o 



-H 



p,„„ = 4p«^^/ {pi + 3), p^''^ = 30o//o78^. 
where pf ^'^ is the critical density of BDA. 



(15) 
(16) 



Using the analogy with the Lemaitre model, Eq. ( flST l is trans- 
formed as 



Q„,„ + Q,A„ + Q^„ — 1 



■^Ao 



'00 



Here, energy density parameters are defined as 

O - P'"o r. _ "^A* - l)Pmo 
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(17) 
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The value Q.^^ is found to be very small as 7.01 x 10"^ for 
fi = 0.6. If we consider the absolute value of O^^ in the parameter 
range 0.0 < n < 0.6, its contribution to Eq. ([TtT i is always less 
than 10"^. Therefore as long as we consider the present epoch, 
contribution from O^^ can be neglected. 

Figure [3] shows the m - z relation in BDK for SNIa obser- 
vations. Matter is dominant in this model. The energy density 
of the cosmological term is always less than 20% in the best 
fit parameter region predicted in Sect. 3. The energy density 
of the cosmological term takes always negative values in the 
obtained parameter region. The parameter B* is not eff'ective 
to change the values of Q.,,^^ and Qa„. Since this model is mat- 
ter dominant, it can not be constrained by the SNIa observations. 

The Friedmann model with the energy density parameters of 
(Q,„,Qa) = (1.0,0.0), is merged with this BDA model having 
reduced ;if^ = x^ 1^ - 4.1 17 {whtxt x^ - 2293 and A^ is defined 
as degrees of freedom). This is inconsistent with the present ac- 
celerating universe, which should contain the suflicient amount 
of dark energy to accelerate the universe. To explain the present 
accelerating universe, it needs some modification to the cosmo- 
logical term. 

As the next approach, BDA is modified by adding another 
constant cosmological term Ac,,. The expansion rate in this 
model is written by 



H 
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4\0 
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The present matter density is 

4(l-A,Jpf^^ 
(// + 3) 



(20) 



(21) 



Here the energy density parameter of the constant cosmological 
term is fixed to be 0.7. 
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Fig. 3. Distance modulus vs. redshift for the flat universe in the 
Friedmann model and the BDA with and without constant cos- 
mological term constrained by SNIa observations (Riess et al. 
120071 Astier et al. i2006l Kessler et al. 12009b . 

We find that this model is consistent with the SNIa observa- 
tions as seen in Fig. [3] Total cosmological term becomes large in 
this model and consistent with the present accelerating universe 
with reduced ;\^2 ^ o.98 (where ;t'^ = 546.92). Foryii = 0.5, BDA 
with A,.„ predicts Oa = -4.285 x 10"^ and Q,„ = 0.343. Qa al- 
ways gets negative values in the parameter region of yu predicted 
in Sect. 3. If we consider the total value of energy densities, con- 
tribution from Qao + ^A, to the total energy density is always 
between 60% - 67%. Therefore the cosmological term is dom- 
inant in the present epoch and it can be constrained from the 
present SNIa observations. It is concluded that BDA with Ac„ 
has nearly the same energy density parameters as the Friedmann 
model with (Q„,, Qa) = (0.3, 0.7). Although the cosmological 
term is not important at the early epoch, it plays very important 
role at the present era. All the parameters inherent in BDA be- 
come insufficient as far as the m - z relation is concerned. 



5. Concluding Remarks 

Previous BBN calculations restricted the parameter range as 
-0.5 < // < 0.8 and -10 < B* < 10 for w = 500 (Nakamura 
et al. |2006j . On the other hand, our large value of w - 10** 
leads to decrease the parameter range of B* (-2 < B* < -2). It 
is oppositely affected the other parameter: 0.0 < i-i < 0.6. The 
models parameters are inefficient in the m - z relations of SNIa. 



It should be noted that from Eq. (|7]), A ~ Pm/<P snd at the 
present epoch, Aq is directly connected with p„g . Dark energy 
is written in terms of dark matter However, dark energy and 
dark matter should be distinguishable to give rise to an accel- 
erated expansion, since evolution of the scale factor seriously 
depends on the composition of each energy density of the uni- 
verse. Therefore, BDA without a constant cosmological term is 
indistinguishable from the matter dominant Friedmann model 
with the parameters of (Q.,„,Q./^) = (1.0,0.0). It is noticed that 
the variable A term in BDA plays a minor role to accelerate 
the universe at the present epoch. Because of this reason, we 
have done a modification to A by adding a constant cosmologi- 
cal term. It has no relation to the expansion rate of the universe at 
the early epoch. However, the energy stored in the constant cos- 
mological term has done a major role to accelerate the universe 
at the present epoch as seen in Fig. [3] Since this model contains 
enough dark energy to accelerate the universe, it is constrained 
by the SNIa observations. In the present research, we have in- 
vestigated BDA at the early epoch to determine the intrinsic pa- 
rameters and introduce new parameters at the present epoch for 
the m - z relation. Since we have demonstrated a possibility of 
non-standard model which is compatible with the observations, 
it is worthwhile to examine more general functional form to the 
cosmological term (eg. Fukui et al l2001l l. 
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In Sect. 4, the value of Q^^ is found to be much smaller 
compared with the other terms in Eq. ( fTTl i. Even though w is 
increased until 10"* the contribution from Q.^„ to Eq. ( [TtI i is 
always less than 1% in the particular parameter range. There 
is no considerable wrong effect from the assumption we made 
in Sect. 4 to neglect the value of Q^^. In the parameter range 
0.0 < n < 0.6, A has taken negative values according to Eq. 
(fTTI ). This may not conflict with theories, since the pressure 
of dark energy must be negative to reproduce the present 
accelerated expansion (Carroli ri992| l. 



